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1.  Introduction 


Anyone  who  is  familiar  with  recent  developments  in  continuum 
mechanics  is  aware  of  the  many  special  features  of  the  equations 
of  hyperelasticity  and  thermoelasticity.  From  the  viewpoint  of 
analysis,  the  most  aggravating  and,  at  the  same  time,  most 
challenging  characteristic  is  that  internal  energy  cannot  be  a 
globally  convex  function  of  deformation  gradient  without  violating  the 
principle  of  material  frame  indifference  [1], 

A  consequence  of  this  lack  of  convexity  is  that  the  classical 
direct  methods  of  the  calculus  of  variations  are  ineffective  for 
establishing  existence,  uniqueness  and  stability  of  equilibrium 
solutions.  Among  the  various  restrictions  on  material  response 
that  have  been  laid  down  for  study  in  the  theory  of  elasticity  [1], 
Strong  Ellipticity  appears  interesting  and  promising.  The  work 
of  Ball  [2]  has  produced  a  breakthrough  by  establishing  the 
existence  of  equilibrium  solutions  in  isotropic  hyperelasticity 
under  the  assumption  that  the  stored-energy  function  is  polyconvex, 
a  condition  somewhat  stronger  than  Strong  Ellipticity  but 
definitely  weaker  than  convexity  and,  in  particular,  not  in¬ 
compatible  with  frame  indifference. 

^^In  dynamics,  the  goal  is  to  relate  the  Second  Law  of  thermo 
dynamics  with  ^stability*.  The  work  of  Ericksen  reveals  that 
the  Second  Law  induces  stability  of  an  equilibrium  state  of  a 
thermoelastic  material  in  the  range  of  convexity  of  internal 
energy  (e.g.  at  points  where  energy  attains  a  minimum.)  Similarly, 
it  has  been  shown  Wrf  that  smooth  adiabatic  processes  in 
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thermoelasticity,  residing  in  the  range  of  convexity  of  internal 
energy,  depend  continuously  upon  the  initial  state,  even  within  the 
class  of  processes  that  develop  shock  waves.  Actually,  in  certain 
circumstances,  the  same  result  can  be  established  under  a  mere 
Strong  Ellipticity  assumption  ^ 

Here  we  investigate  the  implications  of  the  Second  Law  upon 
a  different  manifestation  of  stability,  namely  the  stability  of 
weak  shock  waves  expressed  by  the  E-condition  of  Lax  [5].  If 
energy  were  convex,  the  results  of  Lax  [6]  on  general  hyperbolic 
systems  of  conservation  laws  endowed  with  a  convex  "entropy" 
would  imply  that  the  Second  Law  is  equivalent  to  the  E-condition. 

A  property  of  this  nature  was  established  by  Malek-Madani  [7]  in 
the  context  of  hyperelasticity  under  the  weaker  assumption  of 
Strong  Ellipticity.  In  [7]  an  "energy  criterion"  plays  the  role 
of  the  Second  Law.  The  work  of  Ericksen  [8]  suggests  that  certain 
phase  transitions  may  be  interpreted  through  failure  of  Strong 
Ellipticity  so  it  would  be  advisable  to  avoid  adopting  this  con¬ 
dition  at  the  outset.  Consequently,  our  project  here  is  to  work 
in  the  broader  framework  of  adiabatic  thermoelasticity,  in  which 
the  Second  Law  is  expressed  in  a  definitive  form  by  the  Clausius- 
Duhem  inequality,  and  to  dispense  with  the  need  of  imposing  Strong 
Ellipticity  in  advance.  We  show  that,  due  to  the  special  structure 
of  the  equations  of  elasticity,  the  Second  Law  implies  that  every 
weak  shock  satisfies  Lax's  E-condition,  without  regard  to  convexity 
or  Strong  Ellipticity. 

It  would  be  interesting  to  examine  the  implications  of  lack 
of  convexity  upon  other  shock  admissibility  criteria.  This  program 
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is  carried  out,  to  a  certain  extent,  in  [7],  within  the  context 
of  hyperelasticity,  and  in  [9],  in  the  context  of  one  dimensional 
elasticity. 


2.  Adiabatic  Processes  in  Thermoelasticity 


We  consider  a  body  with  reference  configuration  Si  <=  Rn, 
n  =  2  or  3,  and  reference  density  P(X).  A  thermodynamic  process 
of  Si  is  determined  by  a  motion  x(X,t)  and  a  (specific)  entropy 
field  n(X,t) .  A  motion  generates  a  velocity  field  v  ■  x  and  a 
deformation  gradient  field  F  =  V^x,  det  F  >  0.  An  adiabatic 
process  with  zero  body  force  is  governed  by  the  conservation  laws 
of  momentum  and  energy,  viz. , 


ia 


PE 
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(2.1) 


and  the  Clausius-Duhem  inequality 

pfi  >  0,  (2.2) 

where  T  is  the  Piola-Kirchhoff  stress  and  E  is  the  (specific) 
energy,  sum  of  kinetic  energy  and  internal  energy  €, 
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♦  €. 


(2.3) 


In  thermoelasticity  theory  €,T  and  temperature  9  are 
determined  by  F  and  n  via  constitutive  relations 

€  =  € (F, n;X) ,  (2.4) 

/s 

Tia  =  T.a(F,n;X)  =  P  ,  (2.5) 

A 

0  -  0 (F»n;X)  =  ||  >  0.  (2.6) 

A 

For  every  proper  orthogonal  matrix  R,  the  function  €  should 
satisfy  on  its  domain  the  condition 

€(RF,n;X)  =  e(F,n;X) ,  (2.7) 

which  is  dictated  by  the  principle  of  material  frame  indifference, 
and  possibly  additional  restrictions  expressing  material  symmetry. 

One  of  the  implications  of  (2.7)  is  [1,§52]  that  €  cannot 
be  a  globally  convex  function  of  F,  for  fixed  n  and  X.  Weaker 
than  convexity  and  not  necessarily  incompatible  with  (2.7)  is 
rank-one  convexity  (Hadamard’s  condition).  Uniform  rank-one 
convexity  or  Strong  Ellipticity  requires  that  at  every  point  in 

A  n 

the  domain  of  €  and  for  any  unit  vector  N  in  R  the  matrix 
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is  positive  definite.  This  condition  implies  that  the  system 
(2.1)  is  hyperbolic.  The  characteristic  speeds,  in  the  direction 
N,  are  0,  with  multiplicity  n(n-l)  +  1,  and  ±AT,  i  ■  l,...,n, 
where  are  t*ie  eigenvalues  of  P(N). 

3.  Lax's  E- condition  and  the  Second  Law  of  Thermodynamics 

In  view  of  (2.5)  and  (2.6)  every  smooth  process  that  satisfies 
the  conservation  laws  (2.1)  satisfies  automatically  the  Clausius- 
Duhem  inequality  (2.2),  as  an  equality.  However,  smooth  solutions 
of  (2.1)  generally  break  down  and  shock  waves  develop.  It  is  in 
the  class  of  discontinuous  solutions  that  the  role  of  (2.2) 
becomes  important. 

Functions  of  bounded  variation  in  the  sense  of  Tonelli-Cesari 
[10]  constitute  the  natural  class  in  which  solutions  of  (2.1) 
should  be  sought.  Functions  in  this  class  are  endowed  with 
geometric  structure  that  closely  resembles  that  of  piecewise  smooth 
functions.  In  particular,  within  the  class  of  solutions  of 
bounded  variation  one  may  distinguish  shock  waves  across  which 
the  classical  Rankine-Hugoniot  jump  conditions  are  satisfied, 
namely. 


slFiaJ  +  IviJNa  "  0 

Ps[v.)  *  [Tia]Na  -  0  (3.1) 

PS[E]  +  [Tiavi]Na  -  0 
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where  N  is  a  unit  vector  pointing  in  the  direction  of  propagation 
and  s  is  the  speed  of  propagation  of  the  shock.  In  (3.1)  and 
throughout,  a  bracket  [  ]  denotes  the  jump  (left  minus  right)  of 

the  enclosed  quantity  across  the  shock.  The  jump  condition 
associated  with  (2.2)  reads 


Ps[n]  <  0.  (3.2) 

A  shock  that  satisfies  (3.2)  will  be  called  admissible. 

It  can  be  shown  [11]  that  strict  rank-one  convexity  is 
equivalent  to  nonexistence  of  isentropic  shocks  of  speed  0  (if 
one  allows  [n]  f  0,  this  is  no  longer  true).  On  the  contrary, 
as  the  following  two  propositions  indicate,  positive  eigenvalues 
of  P(N)  and  weak  shocks  of  nonzero  speed  go  hand  in  hand. 

Proposition  3.1.  Assume  that  every  neighborhood  of  a  state 

(F*,v*,n  )  contains  some  point  (F+,y+,n+)  which  can  be 

connected  to  (F',v’,n)  with  a  (not  necessarily  admissible) 

shock  propagating  with  speed  s,  |s|  >  Sq  >  0.  Then  there  is  a 

unit  vector  N  in  Rn  such  that  the  matrix  P(N),  defined  by 

2 

(2.8),  at  (F  ,v  ,n  )  has  at  least  one  eigenvalue  X  >  sQ. 

Proof.  We  fix  a  point  (F+,v+,n+)  which  can  be  connected  to 
(F’,v',n  )  with  a  weak  shock  propagating  in  the  direction,  say, 
N  with  speed  s.  From  (3.1)j,  s2 [Fio]  [Fia]  *  [vjHvj]  so  that 
| [F] |  and  |[v]|  are  of  the  same  order  of  magnitude.  We  now 
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show  that,  in  contrast,  [n]  is  of  higher  order  than  | [v] | . 
Using  (3.1)3,  (2.3),  (2.5)  and  (2.6), 

0  •  PS [EJ  *  lTiaviJNa  ■  J  PstVjVj]  ♦  Ps[£J  *  [Tiav^ ]Na 

-  1  psCv‘.v:)[v1)  .  sT:a[Fla]  .  \  PS  a/fp  [FlaHFjfij 

lu  Jp 


+  PS0"[n]  +  [Tiavi]Nct 


(3.3) 


where 


0-=0(F+,n), 


32€ 


air  am 

ia  j  e 


9a 

"5F!  3F  ~ 
i<*  j  e 


CE,n') 


with  n  between  n  , n  and  F  on  the  straight  line  segment  in 
2  .  +~ 

R  which  joins  F  with  F  .  Hence,  on  account  of  (3.1)  and 

(2.8), 

- p sQ  [ n ]  =  \  ps(v%vT)  [vi]  +  +  lTiavi]^ 

+  ij CS)tvi3tv.l 

*  \  Ps(v%v:)[Vi]  ♦  vt[Tia]Na  ♦  \  f  ^(NJIv.Hv.] 

=  \  ps(v%v:)[vi]  -  psvt[vi]  +  \  f  PijfNJfv.JIVj] 

m~\  pslv.][v.]  +  \  f  Pij(S)[vi][v;.], 


(3.4) 
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which  shows  that  [n]  is  at  most  of  order  | [v] |  . 

On  the  strength  of  the  above  results,  combining  (3.1)^  with 
(3.1)2  we  obtain 

P..(N)[v.]  -  s2[Vi]  =  0(| [v] | 2)  (3.5) 

A 

with  P (N)  evaluated  at  (F’,n).  Dividing  through  by  | [v] | 
in  (3.5)  and  passing  to  the  limit  as  |  [v]  |  -*■  0,  using  the 
compactness  of  the  unit  sphere  in  Rn  we  arrive  at  the  assertion. 

The  proof  is  complete. 

Proposition  3.2.  Assume  that  P(N),  at  a  state  (F°,n°)  and  a 

vector  N,  has  a  simple  positive  eigenvalue  Xq.  Set 

Mo  ■  ^  or  V*0  '  an<*  an^  ~°  Tlien  there  are 

2 

two  smooth  maps  from  an  interval  to  R  and  Rn  +n+1  which  carry 

t  €  (-6,6)  into  s(t)  and  (F,v,n)(t),  respectively,  with  the 

property  s(0)  =  , (F,v, n) (0)  =  (F°,v°,n°)  and  such  that 

(F,v,n)(t)  can  be  connected  to  (F°,y°,ri0)  by  a  shock 

propagating  in  the  direction  N  with  speed  s(^),  i.e., 

s(Fia-F?a)  +  Cvrv?>Na  =  0 

Ps(vrv®)  *  (Ti0-T°0)N0  -  0  (5.6) 

Ps(E-E°)  ♦  (Tiavi'Ti«vi)NP  -  0. 

Furthermore,  (F,v,n)(t),  t  e  (-6,6),  are  the  only  states  in  some 
neighborhood  of  (F0,v°,n0)  that  can  be  connected  to  (F°,y°,n0) 
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by  a  shock  of  speed  close  to  Pg. 


Proof.  First  we  observe  that,  since  P  >  0  and  3(pE)/3n  =  P9  >  0,  the  Jacobian 
matrix  3(F,Pv,PE)/3(F,v, n)  is  nonsingular.  Thus,  we  may  visualize  (3.6)  as  an 
equation  sZ  -  G(Z)  =  0  where  Z  =  (F-F°,Py-Py°,PE-PE0) .  We  observe  that  G(0)  = 
Furthermore,  since  is  a  simple  eigenvalue  of  P(N),  PQ  is  a  single  character¬ 
istic  speed  of  (2.1)  so  that  pQ  is  a  siirple  eigenvalue  of  the  matrix  Gz(0).  Hie 
assertion  of  the  proposition  then  follows  from  a  standard  theorem  in  bifurcation 
theory  [12]. 

We  now  assume  that  the  conditions  of  Proposition  3.2  hold  and 
we  discuss  the  admissibility  of  the  shocks  associated  with  the  shock 
curve  (F,v,n)(T). 

For  t  near  0,  P(N)  ,  evaluated  at  (F,n)(T),  will  have  a 
simple  eigenvalue  A(x)  near  and,  therefore,  (2.1)  will  have 

a  characteristic  speed  p(t)  near  Pg.  The  shock  which  connects 
(F,v,n)(*)  with  (F°,y°,n0)  will  satisfy  Lax's  E-condition  [5]  if 

P(t)  <  s(t)  <  Pq.  (3.7) 

This  condition  may  be  motivated  in  a  variety  of  ways.  For  example, 
it  can  be  interpreted  as  a  stability  statement  for  the  systems 
derived  by  linearizing  (2.1)  on  both  sides  of  the  shock. 

If  (2.1)  were  genuinely  nonlinear  [5]  and  €  were  a  uniformly 
convex  function  of  (F,n)  (in  which  case  -n  would  be  a  uniformly 
convex  function  of  (F,pv,pE)  [13]),  (3.2)  and  (3.7)  would  be 
equivalent  by  a  theorem  of  Lax  [6] .  We  wish  to  investigate  the 
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relationship  between  (3.7)  and  the  Second  Law  of  thermodynamics 

/s 

without  imposing  any  convexity  restrictions  upon  €. 

In  order  to  avoid  the  degenerate  case  where  s'(t)*  changes 
sign  infinitely  many  times  in  every  neighborhood  of  0,  we  make  the 
assumption 


s ^  J(0)  f  0,  for  some  i  >  1 , 


(3.8) 


which  is  generically  satisfied.  For  8,  =  1,  (3.8)  expresses  the 
genuine  nonlinearity  of  the  Pq  characteristic  field  at 
(F°,y0,n°)  and  it  will  be  satisfied  if  and  only  if 


~3~,  co  o. 
3  €  (F  ,  n  ) 


r,  t  0, 


3Fia3Fj63FkY  «  Pt  i  j  k 


(3.9) 


where  r  is  the  eigenvector  of  P(N)  associated  with  the  eigenvalue 


We  are  now  prepared  to  state  the  main  result: 


Proposition  3.3.  Under  assumption  (3.8)  and  for  *  sufficiently 
close  to  0,  Lax's  E-condition  (3.7)  is  equivalent  to  the  strict 
Clausius-Duhem  inequality 


s(*)(n(*)  -  nu)  <  0. 


Here  and  throughout  a  prime  (')  denotes  derivative  with 
respect  to  x. 


(3.10) 
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Proof.  Differentiating  (.3.6)  with  respect  to  x, 


sF! 

+ 

v.’N^  = 

ia 

1  ci 

Psv! 

+ 

T !  N  = 

i 

ia  a 

psE’ 

+ 

T!  v-  N 
ia  i  a 

-s ' (F.  -F:  ) 

L  iu  iaJ 


Using  (2.3),  (2.5)  and  (2.6),  (3.13)  yields 


psvivi  +  sTiaFia  +  ps0r1'  +  TiaviNa  +  TiaviNa  =  -ps'(E-E°; 


whence,  by  virtue  of  (3.11)  and  (3.12), 


PsOn’  =  -s ’ { P (E-E°)  -  Tia(Fia-F°a)  -  Pvi(vi-v°)) 


On  account  of  (3.6),  (3.14)  gives 

PsW  -  <■  Tia(vrv°)N0  -  psvjtvj-v”) 

=  -s''-'ri(Tia-Tia’N<»  ‘  PsVvry°» 

=  -s’ (PSV^Vj-V?)  -  Psvi(vi-v?)} 


=  Ps  s 1 (Vj-v?) (vi-v?) 


(3.11) 

(3.12) 

(3.13) 


(3.14) 


Hence , 
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sen'  =  S' (vi-v?)(vi-v°).  (3.15) 

In  view  of  (3.8),  for  x  near  0,  s'  does  not  change  sign 
between  0  and  x.  it  then  follows  from  (3.15)  that  n'  also 
does  not  change  sign  between  0  and  x  and  that  sn’  and  s' 
have  the  same  sign  on  this  interval.  Thus,  (3.10)  is  equivalent 
to  sCx)  <500)  =  Pq,  X  e  (0,T).  By  reversing  the  roles  of  left  and  right 
state  one  shows  that  (3.10)  is  also  equivalent  to  s(x)  >  y(x).  The  proof 
is  complete. 
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